Spatial effects in superradiant Rayleigh scattering from Bose-Einstein condensates 
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We present a detailed theoretical analysis of superradiant Rayleigh scattering from atomic Bose- 
Einstein condensates. A thorough investigation of the spatially resolved time-evolution of optical 
and matter-wave fields is performed in the framework of the semiclassical Maxwell-Schrodinger 
equations. Our theory is not only able to explain many of the known experimental observations, e.g., 
the behavior of the atomic side-mode distributions, but also provides further detailed insights into 
the coupled dynamics of optical and matter-wave fields. To work out the significance of propagation 
effects, we compare our results to other theoretical models in which these efi'ects are neglected. 
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I. INTRODUCTION 

The beautiful recent experiments of Refs. P, Q have 
opened up the possibility of studying superradiance 
within the context of Bose-Einstein condensates of ul- 
tracold atomic gases. In a typical experimental setup, 
a cigar-shaped condensate is exposed to a far-off reso- 
nant laser pulse travelling in a direction perpendicular 
to the long condensate axis [see Fig. ^a)] . Condensate 
atoms can then undergo Rayleigh scattering thereby ex- 
periencing a recoil kick. The moving atoms together with 
the condensate at rest form matter-wave gratings from 
which further laser photons are scattered. This in turn 
produces additional recoiling atoms, causing the ampli- 
tudes of the matter- wave gratings and the scattered light 
fields to grow rapidly in a self-amplifying process. Due 
to phase-matching effects induced by the elongated con- 
densate shape, the fastest growth is experienced by those 
gratings for which the scattered photons travel parallel 
to the condensate axis in the two so-called optical "end- 
fire modes." After an initial start-up and mode competi- 
tion process 0, the ensuing superradiant light emission 
is concentrated into the endfire modes, and the recoil- 
ing atoms have well-defined momenta thus forming two 
first-order atomic "side modes" [Fig. Q^a)] . 

Comparing this process to "conventional" superradi- 
ance in an inverted atomic medium 0, 1^, we see that 
the role of the electronically excited atoms is now taken 
by the combination of the condensate and the applied 
laser pulse, spontaneous photon emission is replaced by 
Rayleigh scattering, and de-excited atoms correspond to 
atoms in the momentum side modes . Another impor- 
tant difference to the conventional process, which stops 
after each atom has undergone a single transition, is the 
possibility of repeated or modified scattering cycles which 
lead to the population of multiple atomic side modes. 
These processes occur naturally if the applied laser pulse 
is of sufficient duration or strength. On the one hand, 
an atom in a side mode can again scatter a laser photon, 
thereby transferring to a higher-order side mode. Alter- 
natively, however, atoms may also interact with endfire 
mode photons and scatter them back into the laser field 
0. This leads to the production of atoms moving back- 



wards with respect to the direction of the applied laser 
field. Since the latter process requires to overcome an 
energy mismatch, it only occurs for sufRciently intense 
laser pulses. In this way, one is led to distinguish be- 
tween the weak- and the strong-pulse regimes of super- 
radiance. The behavior of the system strongly differs in 
these two regimes. In particular, one important differ- 
ence concerns the atomic side-mode population patterns 
which can be investigated through time-of-flight imag- 
ing. In the weak-pulse regime, where forward scattering 
is prevalent, distributions take a characteristic fan shape, 
whereas in the strong-pulse regime, an X-like pattern is 
observed. 

Basic features of the atomic momentum distributions 
were theoretically discussed by a number of authors (see, 
e-g-, [E IE Hi E) without, however, giving a com- 
prehensive description of the experimental observations. 
Reference [Tol |. for example, examined a model in which 
two regimes of either pure forward scattering or combined 
backward and forward scattering could be distinguished 
depending on the external parameters. However, this 
study only considered a strictly one-dimensional system 
so that a direct comparison to the experiments was not 
possible. 

A more comprehensive investigation of the atomic side- 
mode distributions was recently presented in our work 
[Tl] |. In Ref. 11], the typical X- and fan-shape patterns 
along with some of their characteristic properties were re- 
produced and explained in terms of the underlying cou- 
pled dynamics of optical and matter- wave fields. Our 
approach is built upon two main concepts: First of all, 
we use a semiclassical description for the field dynamics. 
Follo wing earlier studies of conventional superradiance 
[E0,0|, we can expect that the semiclassical approach 
is valid as soon as the numbers of atoms and photons 
in the atomic side and optical endfire modes, respec- 
tively, become large compared to one. It is this macro- 
scopic regime of superradiance that we are focussing on 
throughout this work. A fully quantized model, on the 
other hand, such as the one used by Meystre and co- 
workers U Q, allows to investigate the initial startup 
of the superradiant process, but is not easily extended 
to the study of long-time dynamics. The second main 
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feature of our approach is the inclusion of spatial propa- 
gation effects. This aspect has been neglected in almost 
all previous quantal as well as semiclassical treatments 
of Bose-Einstein condensate (BEC) superradiance (see, 
e.g., ,3, 6j, 2s, ^ Me JM)- References Il5| examine 
some spatial effects in the BEC-light interaction, but do 
not provide a detailed comparison to the experimental 
observations of Refs. [Hi]. The results of Ref. Tl], how- 
ever, show that the inclusion of spatial effects is crucial 
for a full understanding of BEC superradiance. In fact, 
our results also allowed us to resolve the controversy be- 
tween Refs. and regarding the origin of the spa- 
tial asymmetry between forward- and backward-moving 
atomic side modes observed in the strong-pulse regime of 
superradiance. 

This successful explanation of several key observations 
of the superradiance experiments 0, Q suggests that the 
model of Ref. ^T^\ provides an adequate description of 
the system dynamics in the semiclassical regime. The 
purpose of the present paper is therefore to extend our 
previous work by providing an in-depth examination of 
this model and, in this way, to obtain further detailed 
insights into the system dynamics. We also continue to 
stress the significance of spatial propagation effects by 
explicitly comparing our results to a related mean-field 
model in which these effects are neglected (see, e.g., [Toj). 

The paper is organized as follows. In Sec. II, we present 
our theoretical framework which is based on the semi- 
classical Maxwell-Schrodinger equations and the use of 
the slowly-varying-envelope approximation. We also in- 
dicate how one may include photon scattering within or 
between the optical endfire modes. Section III analyti- 
cally discusses the short-time gain regime. In particular, 
it is shown that due to the propagation effects, the growth 
of the atomic side-mode populations and optical endfire- 
mode intensities is slower than exponential. In Sec. IV, 
we discuss the strong-pulse regime and explain in de- 
tail how the characteristic observations of X-shape side- 
mode patterns and the spatial asymmetry come about. 
We also exhibit explicitly the failure of spatially indepen- 
dent models to account for these effects. In Sec. V, the 
weak-pulse regime is investigated. Besides discussing the 
atomic side-mode distribution patterns, we also use our 
model to study characteristic spatial effects in the cou- 
pled dynamics of optical and matter- wave fields. It will 
turn out that these effects are significantly more involved 
than in the strong-pulse regime. The paper concludes 
with a short summary and outlook in Sec. VI. 



II. MAXWELL-SCHRODINGER EQUATIONS 

Our model involves an elongated condensate of length 
L oriented along the z axis. The BEC consists of 
N two-level atoms which are coupled via the electric- 
dipole interaction to a linearly polarized pump laser pulse 
£i(i)ey (e^^*^'^-"'*) -t-c.c.)/2 with = c/c;, traveling in the 
X direction [see Fig. BJa)]. The laser is considered to be 




FIG. 1: Schematic representation of BEC superradiance. (a) 
A cigar-shaped BEC of length L (filled ellipsoid) is exposed to 
a linearly polarized laser pulse with amplitude £i (t) and wave 
vector ki. Stimulated Rayleigh scattering of laser photons 
off the condensate produces recoiling atoms in side modes 
with well-defined momenta. The side mode (0, 0) refers to 
the BEC at rest whereas only one of the two first-order side 
modes, i.e., (1, —1), is shown. The scattered photons propa- 
gate mainly along the condensate axis in the endfire modes 
£±. (b) Population transfer to and from the side mode (n, m) 
due to interaction with the laser pulse as implied by Eq. (|7|l. 



far off resonant from the atomic transition to the excited 
electronic state |e). 

As discussed in the Introduction, in this system the 
coherent nature of the BEC leads to strong correlations 
between successive Rayleigh scattering events and to col- 
lective superradiant behaviour Q,0,y{. Moreover, as a 
result of the cigar shape of the condensate, the gain is 
largest when the scattered photons leave the condensate 
travelling up and down its long axis in the so called end- 
fire modes £±(x, f)eye~**^"*^'^^\ where £±(x, t) are the 
envelope functions of the modes. As a consequence, the 
recoiling atoms have well-defined momenta and appear 
in distinct atomic side modes. In the side mode (n,m), 
atoms possess momentum h{nkiex + mkez), while their 
kinetic energy is given by hujn,m = ft^ {n'^ kf + m'^ k^) / 2M . 
In this notation, the "side mode" (0, 0) describes the con- 
densate at rest. The wave vectors ifce^ of the scattered 
photons are fixed by energy conservation for the tran- 
sitions between the side modes (0,0) and (1,±1) which 
initiate the process, i.e., hcki = hck + hwn. Given that 
^n,m ijJi, we approximately have k ^ ki and ki — k <^ 
k,ki. Thus, the side-mode frequency is approximately 
given by ujri,m ~ {n^ +m'^)ujr, where tOr = hkf /2AI is the 
recoil frequency. 

Maxwell-Schrodinger equations of motion for the cou- 
pled matter-wave and electromagnetic fields offer a very 
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general theoretical framework for the description of prob- 
lems per tainin g to the interaction of ultracold atoms with 
light ■f'or the problem at hand, after adiabati- 

cally eliminating the excited atomic state |e), the coupled 
Maxwell-Schrodinger equations for the mean- field macro- 
scopic atomic wave function ')p{'x., t) and the positive- and 
negative-frequency components E(='=)(x, t) of the classical 
electric field read 5, 16, 17] 



in—w = Alb 

dt ^ 2M ^ 

52e(±) 



(d-E(-))(d •£(+)) 



9t2 



1 6>2P(^) 
9t2 



^, (1) 
(2) 



with 5 the detuning of the laser from the electronic tran- 
sition, d the atomic dipole moment and M the atomic 
mass. The polarization is given by 



p(+)(x,i) = -d|V'(x,i)| 



2 d-E(+)(x,t) 

h5 



(3) 



with P(") = p(+)*. Note that in Eq. (jT} we neglect 
the external trapping potential and interactions between 
atoms, since they do not play a significant role on the 
time scales of the process under consideration. However, 
it would be straightforward to include them in the model. 

In order to solve Eqs. ^ and l(2Jl under the slowly- 
varying-envelope approximation (SVEA), we decompose 
the fields as 

V'(x,i) = ^"'"^'^^ e-'^""^"'*-"^'"-"'''-), (4) 

/ N w A. 
[71, m) 

E(+)(x,i) = £ieye-'^'^''~'"''^/2 + £+{z,t)eye-''^'^'-'"'> 
+£-{z,t)eye~'^'^'+'''\ (5) 



E(~)(x, t) = E(+) *(x, t), with u! = kc and A the average 
condensate cross section perpendicular to the z axis. The 
electric field in Eq. © contains the impinging laser pulse 
together with the two optical endfire modes that are pro- 
duced by collective Rayleigh scattering. The atoms in 
the side mode (n, to) are characterized by a slowly vary- 
ing spatial envelope ipnmiz, t). The summation in Eq. 
can be restricted to terms (n, to) with m-\-n even, as Eq. 
(UI) below will show. 

Note that in the ansatz Q and (0) we disregard the 
dependence of the envelope functions ^nm and £± on the 
transverse directions x and y. For the matter waves, this 
is certainly a good approximation since the radial degrees 
of freedom are tightly confined by the trap. Moreover, 
the use of the same approximation for the optical fields is 
justified by the fact that in the ex]Deriments the Fresnel 
number of the system is close to 1 
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In order to introduce a concise notation, we rescale the 
optical and matter waves as 



hhjki 
2^' 



Ipnr. 



A 



(6) 



Then, using the ansatz Q and lO and introducing the 
dimcnsionless time t — 2ujrt and length ^ = fc/z, Eq. 
in the SVEA reads 



e; (C, r) V'„_i,™+i (e, r)e*("-'"-2)- + e*_ (C, t)^„_i.™_i (C, r)e'("+'"-2)- 
r)V'„+i,™-i(e, r)e-^("-™)^ + e_(e, T)V'„+i.„+i(e, r)e-^("+")^' 



-A 



el (e, r)e+(e, T)^„,„_2(e, T)e^^^"^-^^^ + e;(?, T)e_ (C, T)^„,„+2(e, r)e-2.(™+i)r 



-A(|e+(e,r)|2 + |e_(e,r)|2)V'„„(e,r), 



r 



(7) 



with the coupling constants 



g 

2uJr 
K 

eo' 



(8b) 



where 



9 ■ 



\d?£i 



2h'^S V 2eaAL' 



(9) 



The first term on the right-hand side of Eq. Q) describes 
the quantum-mechanical dispersion of the envelope func- 
tion, while the second one leads to a spatial translation 
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with velocity Vm = mhk/M. The other terms describe 
the interaction between the matter- wave and electromag- 
netic fields. This interaction leads to spatially dependent 
shifts and couplings of the momentum side mode (71, m) 
to other modes. Let us discuss the underlying physical 
processes in more detail. 

The terms involving the coupling constant k refer to 
photon exchange between one of the endfire modes and 
the laser beam [see Fig. H^b)]. In particular, through 
stimulated scattering, an atom in a side mode (n, m) can 
absorb a laser photon and deposit it into one of the end- 
fire modes. The accompanying recoil transfers the atom 
into one of the side modes {n + 1, m ± 1). Alternatively, 
the atom may absorb an endfire-mode photon and emit it 
into the laser beam, thereby ending up in the side mode 
(n—1, TO±1). This latter process is responsible for atomic 



backward scattering 

Another class of physical processes involves the photon 
exchange between the two side modes and is described by 
the terms with e±e^. In this case the accompanying re- 
coil transfers the atom from (n, to) to the side modes 
(n, m±2). Finally, the terms containing |e±(^,T)p de- 
scribe the absorption of a photon from an endfire mode 
and its subsequent emission into the same mode. Hence, 
such processes give only rise to shifts which depend on 
both time and space, but do not introduce couplings to 
other side modes. The corresponding term for the laser 
beam is related to a constant ac Stark shift and as such 
is not included in Eq. (0. 

In the SVEA, the envelope functions e± of the endfire 
modes obey the equations 



'd7 



de_ 



-Z ^ p^"-'")-Vn™(C,T)C+l,„.-l(^,T) 
(n,m) 



(n,m) 



-|-Ae+(e,r)e2*(™+i)-^„„(e,r)V'* +2(?,T)+Ae_(e,T)|V'„™(C,T)|2 



(10) 



(11) 



where 



X 



cki 



(12) 



Introducing the atomic natural linewidth = 
d^a;^/(37r£o^c'^), the couphng constant g can be ex- 
pressed in terms of the Rayleigh scattering rate R — 
(PEfVaKAM'^) as 



37rc3i? 



(13) 



Moreover, the collective superradiant gain is given by 



G 



7 



(14) 



with the photon damping rate 7 = 2c/L (note slight dif- 
ferences from Rcf. 10] regarding the numerical prefactors 
in the definitions of G and 7) . Experimental observations 
clearly distinguish between two different regimes of pa- 
rameters, the so-called strong- and weak-pulse regimes 
P, 0, ■ The former regime is characterized by a scatter- 
ing rate comparable to the recoil frequency and G ^ ujr, 
while for the latter R<^ LOr and G < Ur- 

For both of these regimes, one can verify that photon 
exchange between endfire modes can be disregarded in 



our model (at least for the time-scales of interest) P, | 
Hence, Eq. Q may be simplified to 



n 



(15a) 



dr 2 



-e+^pn 



Since the interaction between the optical and the matter- 
wave fields is restricted to the condensate volume, any 
relevant retardations are expected to be of the order of 
L/c ~ 10""'^^ s and can be neglected. Thus, formal inte- 
gration of Eqs. H10|I - H11|) for A = yields the following 
equations for the envelope functions e±: 



e+(e,T) = 



e-(e,T) 



K 

—i — 
X 



(n,m) 



xV'«m(^', T)V'*,+ l,m-l(^'' (15b) 



K 

~i— 
XJi 



.m 
(n.m) 



xV'„™(r,r)C+i^„.+i(^',r). (15c) 

From these equations one may clearly see how the 
buildup of the endfire-mode fields at (^, t) is driven by 
the coherences V'nmV'^+i,m±i- 
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As mentioned in the Introduction, previous theoretical 
studies of BEC superradiance have neglected the spa- 
tial dependence of optical and matter- wave fields. It is 
one of the main purposes of this paper to compare the 
model derived above to this simpler approach in order 
to assess the significance of propagation effects. For- 
mally, we can obtain the equations of motion for the 
spatially independent model from Eqs. (|15|l by dropping 
any spatial dependence and setting ipnm — > Cnm/ VkiL 
and e± — *■ b±/^/kiL. The amplitudes Cnm{T) and 6±(t) 
can be interpreted as amplitudes of the matter- wave and 
optical fields, respectively, and they obey the equations 
of motion 





, 9 


dr 












7 


n,m) 


&_(r) = -i^ 

7 


E' 

n.m) 



blCn- 



■i,m+ie 



i(7i — m — 2)T 



t(n—m It 



(16a) 

e*("-'")^C„™(T)Q+i,^_i(r), (16b) 
e^("+")^C„™(T)C:+i,„,+i(T). (16c) 



A similar model, which distinguished only the momen- 
tum side modes in the x diretion, was presented in |in| |. 
It should be emphasized that one cannot consider C„m (t) 
and 6±(r) as spatial averages of ipnmi£,,T) and e±(^,r). 
In fact, we will see that the predictions of the two models 
show significant differences. It should also be noted that 
the damping coefficient 7, although of the order of c/L, 
cannot unambiguously be ascribed a precise value. 

In general, superradiant scattering is initiated by 
quantum- mechanical noise, i.e., spontaneous Rayleigh 
scattering from individual condensate atoms Sub- 
sequent stimulated scattering and bosonic enhancement 
lead to rapid growth of the side-mode populations. The 
semiclassical model derived in this section describes the 
macroscopic stage of the superradiant process where the 
populations of the side modes are already large compared 
to one. 

Adapting the discussion of Refs. |^ 0, 0| regarding 
"conventional" superradiance, we can take the effects of 
the initial quantum-mechanical fluctuations into account 
by solving the semiclassical equations of motion with 
stochastic initial conditions (seeds) for the side modes 
V'i,±i(Ci''' = 0)- Since the noise is practically relevant 
for the initial population of those modes only, we can set 
'0nm(Ci T = 0) = for all other side modes. The conden- 
sate '000 is chosen to be in the macroscopic ground state 
at t = 0. For any of these stochastic initial conditions, the 
solution of the semiclassical equations of motion corre- 
sponds to one possible realization of the experiment. By 
studying a large set of simulations with varying seed func- 
tions drawn from an appropriate distribution, one could 
obtain information on, e.g., averages and fluctuations of 
relevant experimental observables. In the present paper, 
however, we intend to focus on characteristic features of 



the macroscopic dynamics that are observed in individ- 
ual experimental realizations, such as those reported in 
Refs. 1113. 

To this end, we have numerically solved Eqs. (|15|l with 
fixed external parameters for a large set of different seed 
functions, and we find the characteristic qualitative fea- 
tures, e.g., the side-mode distribution patterns, to show 
up independent of the choice of the specific initial con- 
dition. This shows that these features are caused by the 
macroscopic dynamics and are not related to quantum 
fluctuation effects. In Sec. IIIII we will explain how this 
inscnsitivity to the initial conditions comes about. 



III. EARLY STAGE OF SUPERRADIANT 
SCATTERING 

In this section, we investigate the early stages of the 
superradiant process in detail. More precisely, we wish 
to examine the regime where the populations of the first- 
order side modes still remain far below the number of 
atoms in the condensate, but are large compared to 
1, so that the semiclassical model is applicable. The 
undepleted-pump approximation for the condensate can 
then be invoked, which allows to derive analytical results 
for the spatial distributions of the side-modes and their 
total populations. With these results, we can (i) compare 
the growth of the side modes in the weak- and strong- 
pulse regimes, (ii) work out differences to the spatially- 
independent model introduced at the end of Sec. II, and 
(iii) explain why the time-evolution of the side modes is 
rather insensitive to the details of the initial seed func- 
tion. The most interesting result of these studies is prob- 
ably the observation that the propagation effects lead to 
a subexponential growth of the side-mode populations 
which is in contrast to the spatially independent model 
where an exponential growth is predicted. 

In the startup regime described above, only the first- 
order side modes for forward- and backward-scattering 
become populated significantly, so that we can restrict 
Eqs. H15|l to these modes and the condensate. Fur- 
thermore, the depletion of the condensate is assumed 
to be neghgible so that wc set ipnoiS.,''') ~ ipooi^jO) 
(undepleted-pump approximation). The equations for 
the sets of modes (±1,±1) and (±1,^1) then decouple. 
Neglecting the effects of free propagation, which is well 
justified for short times, the time evolution of the side 
modes (±1,=f1), for example, is governed by 



« = -000(0 / [e 0_i_i(^,t) 

X J —00 



dr 



xV'oo(e') + ^o*o(e')V'i,-i(e',T)], (17a) 



« = 000(0 / rf^ |.0-ia(C ,t) 

X J —00 



dT 



x0o*o(C')+e''^0oo(e')V'iVi(^''ttb) 
Before discussing the solutions of Eqs. (|17|l . let us briefiy 



6 



consider the corresponding startup regime in the spa- 
tially independent model. Setting Coo(t) = \/N and 
C_i4(r = 0) = 0, one obtains from Eqs. lfTH|l 

4(2 - Ai) 

(18a) 

C-rAr) = ^il^r (e^- - e^-) (18b) 



with C_i.i(r) 
rates Ai = i—^/^^ 



e^^^Cli i(t), Y = G/uJr, and the growth 
iF, A2 — — iT. In the weak- 

pulse regime of F ^ 1, one has Ci,_i(t) « Ci,_i(0)e'"'^/^ 
and |C_i4/Ci,_ip ~ F^/16. For strong pulses with 
F 3> 1, on the other hand, one finds Ci__i(t) w 
Ci,_i(0)v/r7T6 e(i-')Vr75r+W4 and - 1 

(see also 10]). We thus see that the spatially indepen- 
dent model always predicts an exponential increase of 
the side- mode populations, but with growth rates vary- 
ing between F and VT. 

To facilitate the comparison with these predictions, we 
examine the spatially dependent equations (|17|) for con- 
stant initial conditions V'l, -1(^70) = '00, 0-1,1(^,0) = 
and a homogeneous condensate V'oo(C) = \/N/A, < 
f < A. Here, A — kiL denotes the scaled condensate 
length. As shown in the Appendix, Eqs. (|17|l can then 
be solved approximately by Laplace transform and sub- 
sequent inversion through the saddle-point method. In 
the weak-pulse limit, one obtains 



■01 -i(C,t) 



: exp 



so that the side-mode population grows like 



(2v/lV^) , (19) 



7Vi,_i(t)= / df|0l^„l(^,T)| 



^^^^-^(«)expf4VlV 



SttFt 



with the initial population iVi^_i(0) — |0opA. 
In the strong-pulse regime, wc find 



(20) 



0i,-i(e,T) 



■006 



57ri/12 



re 



1/6 



1 



22/3^6^ \ AJ t2/3 

X exp (it + 3e-*^/'^(Fr2^/2A)i/^) .(21) 



The total side-mode population is approximately given 

by 



AT ( ^ ^1,-1(0) 1 



pl/3^2/3 



(22) 



In contrast to Eqs. (|18|) . equations (|20|1 and (|22|l show 
that propagation effects lead to a subexponential growth 
in the side- mode populations: for weak pulses, iVi__i in- 
creases as exp(y'r)/r, whereas for strong pulses, it grows 
like exp(r2/3)/T2. Moreover, from Eqs. (fT^ and (PT|l . 
we see that the side modes also grow subexponentially in 



space, but again display different behaviors in the weak- 
and strong-pulse limits. 

For a spatially varying condensate wave function, e.g., 
a BEC ground state in a trap, it becomes more difficult 
to derive analytical results for the side modes. However, 
the numerical solution of Eqs. (|17|l with the condensate 
having the Thomas-Fermi shape yields side-mode popu- 
lations that are very close to the ones for a homogeneous 
condensate with the same atom number. We thus expect 
laws similar to Eqs. (|20|l and (|22|l to apply in this case 
as well. 

The study of Eqs. H17|l also allows us to understand 
why, after an initial transient, the spatial shape of the 
side-mode wave functions becomes rather insensitive to 
the details of the initial seed function. To this end, we 
consider the Fourier decomposition of an arbitrary initial 
seed 



0'i,-i(e,o) 



c„ exp(i27rne/A). 



(23) 



Since Eqs. ((T7|l are Hnear, 0i^_i(e,T) is obtained 
as the linear superposition X^n '^nV'i'-i('?j where 

■0i"li(C,O) = exp(i27rne/A). Considering the weak-pulse 
regime for concreteness, the Laplace transform method 
shows that, after some time, the wave functions V"!™-! 
behave like 



exp 



(2yiVe7A' 



27m 



47r( yiVe/A-z— e/yiVe/A 



From this result, we can draw two conclusions, (i) All 

wave functions ■0i"lx acquire more or less the same shape, 
i.e., their moduli grow subexponentially in space, whereas 
their phases vary only slowly. We can thus expect their 
superposition to qualitatively show the same behavior, 
(ii) The contributions of higher-order Fourier components 
are suppressed by a factor of approximately. In fact, 
for growing |n|, we find numerically that it takes increas- 
ingly more time to reach the form H24() . and before this 
time is reached, the suppression is even stronger. These 
observations explain the insensitivity of the side-mode 
wave functions to the details of the initial seed. 



IV. STRONG-PULSE REGIME 

In Sees. IV and V, we will use the theoretical model 
developed in Sec. II to present a thorough discussion 
of the coupled dynamics of optical and matter waves 
in the strong- and weak-pulse regimes of superradiance. 
For the numerical calculations, we focus on the exper- 
imental data of Refs. P, In particular, we con- 
sider a s'^Rb BEC of iV = 2 x 10^ atoms, with length 
L = 200/im and cross-section diameter 15/im. The 
condensate is in the Thomas-Fermi regime, so that we 
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can model its wave function as ipooiz) = \/n{z) with 
n{z)=C[{L/2f-z^]Q{L/2-\z\),C = ?,N/AL^. As dis- 
cussed in Sees. II and III, the results are not significantly 
influenced by the shape of the seed function. Hence, for 
the sake of simplicity, throughout our simulations the 
first-order atomic side modes are seeded according to 
'0i,±i(2:,O) = 'iIjoo{z) / y/N , which corresponds to one de- 
localized atom in each of the modes. The applied laser 
pulse is modeled as rectangular lasting from t = up to 
t = tf. Equations H15|l are then propagated using a split- 
step algorithm 19]. The grid {n,m) for the side- mode 
orders is chosen sufficiently large, so that the population 
of the highest-order side modes remains negligible at all 
times. The chosen values for the external parameters, 
such as pulse strength and duration, are comparable to 
those used in the experiments 1, 2]. 

Let us first consider the regime of strong laser pulses 
which is characterized by Rayleigh scattering rates R 
comparable to LUr- In this regime of parameters, ex- 
periments with short laser pulses have shown forward- 
as well backward-recoiling atoms, forming a very char- 
acteristic X-shape pattern. Moreover, there is a notice- 
able spatial asymmetry between forward and backward 
peaks. According to Ref. 0, these observations suggest 
a purely optical picture of superradiance in which atoms 
are diffracted from the optical grating formed by the end- 
fire and pump modes. The intensity of the endfire modes, 
and thus the optical grating, are assumed to be peaked at 
the edges of the condensate. Our model has enabled us 
to reproduce the experimental observations and to verify 
this conjecture. 



A. Side-mode patterns 

In Fig.|2Ia), we display a snapshot of the atomic spatial 
distribution after applying a strong laser pulse to the 
condensate followed by a free propagation for a time tp 3> 
tf. Since we work with a one-dimensional model, we 
calculate the displacement Ax between the condensate 
and the first-order side modes in the x direction as Ax = 
Vrtp with the recoil velocity Vr = hk/M = 5.9 x 10""^ m/s. 
Our result clearly reproduces the asymmetry observed in 
Fig. 2 of Ref. 0]. 

The physical mechanism behind this asymmetry can 
be understood from the spatially resolved dynamics of 
the optical and matter- wave fields as described by Eqs. 
(|15|l . For the sake of illustration, the spatial distribution 
of the condensate, the first-order atomic side-modes and 
the optical-field modes at the end of the strong pulse 
(i.e., at time tf) are plotted in Fig. I^Jb). Clearly, the 
atomic side-modes and the optical-field modes are well 
localized near the condensate edges. This preferential 
growth of the optical fields can be explained in the con- 
text of Eqs. Ijl5b|l and H15c|l which imply that (at least 
for short times) the electric fields £±{z, t) grow monoton- 
ically in the z and — z directions, respectively, and are 
strongest at the ends of the condensate. As a result, ac- 
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FIG. 2: Spatially dependent model, (a) Spatial distribution of 
the first-order forward (1, ±1) and backward ( — 1, ±1) atomic 
side modes, after applying a laser pulse of duration tf = 8.5 
US and strength g = 3 x 10^ s~^ to the condensate followed 
by a free propagation for a time tp — 25 ms. (b) Spatial 
distributions of the atomic side modes and the optical endfire 
modes £± at time tf. For the sake of illustration, in both 
cases the population of the BEG (0, 0) has been multiplied by 
0.25. 



cording to Eq. (|15a|l . the scattering process is strongest in 
these areas of large electric fields, and the recoiling atoms 
mainly originate from the edges of the condensate. Dur- 
ing the free time evolution following the strong pulse, 
the forward-scattered atoms in the side modes (1,±1) 
initially travel towards the center of the BEG, whereas 
the atoms in the backwards side modes (— 1,±1) imme- 
diately move away from the center [dashed lines in Fig. 
|2Ia)]. The net effect is the observed asymmetry in the 
spatial distribution of forward and backward peaks. So, 
in agreement with the experiments, our model predicts 
enhanced photon scattering near the edges of the con- 
densate in the strong-pulse regime. This effect strongly 
supports an interpretation of the strong-pulse regime in 
the framework of atomic scattering from optical fields as 
was first conjectured by Ketterle and co-workers 0. 

Our model also explains the experimentally observed 
X-shape patterns [see Fig. 1(A) of Ref. 2]]. In Fig. 
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and this is possible only if the atomic side modes and 
the optical endfire modes are located at the condensate 
edges. To make this point clear, let us consider, for ex- 
ample, the ofF-diagonal side mode (2,0). Figure 0] shows 
the possible population transfers among the low-order 
atomic side-modes and the optical fields. According to 
Eq. H15a|l . although the side mode (2,0) is resonant with 
the modes (1, ±1), it can only be populated if the (1, ±1) 
modes overlap with the endfire modes £±, respectively. 
As depicted in Fig. El^b), however, this overlap is very 
small since the modes are localized near the edges of 
the condensate, and the growth of the side mode (2, 0) 
is therefore suppressed. On the other hand, population 
transfer to the side modes (2, ±2) is easily accomplished 
due to the strong overlap between the side modes (1, ±1) 
and the fields respectively. 



FIG. 3: Spatially dependent model. Atomic side-mode dis- 
tributions in the strong-pulse regime for g = 3.0 x 10® s~^ 
and tf — 12 us. The gray level of each square represents the 
"relative" probability p„m as defined in the text. The arrow 
indicates the direction of the incoming laser pulse and points 
towards the BEG "side mode" (0,0). 



• (2,2) 




(-1^-1) • (1,-1) • 



• (2,-2) 

FIG. 4: (Color online) Origin of the X-pattern: schematic 
representation. Arrows show population transfer between dif- 
ferent side modes (n, m) via couplings to the optical fields £±. 
The side mode (2, 0) does not become populated due to the 
missing spatial overlap between the side modes (1,±1) and 
the endfire modes £±, respectively. 



13 we show the atomic side-mode distribution forming 
an X-shape pattern at the end of a strong pulse with 
g = 3 X 10^ s~^. The gray level of each square repre- 
sents the "relative" probability Pnm = Pnm/ Pmn^'^ with 
/ dz\'ii)nm{z,t)\'^ and As 
a matter of fact, the enhanced photon scattering near the 
edges of the condensate is also responsible for the charac- 
teristic X-shape distribution of the side modes. Indeed, 
the appearance of an X-shape pattern requires the sup- 
pression of all off-diagonal side modes with \n\ ^ \m\ 



B. Discussion 

Our studies show that the experimentally observed 
asymmetry and X-shape patterns are basically due to the 
spatial inhomogeneity of the optical fields. In an earlier 
attempt to theoretically describe these observations, Pu, 
Zhang, and Meystre (3] adopted a model which does not 
include any spatial propagation effects. In the framework 
of this model, they concluded that the observed asymme- 
try is due to angular effects, thus questioning the picture 
of atomic scattering from optical fields suggested by Ket- 
terle and co-workers 0. In particular, they found that 
forward recoiling atoms should have a symmetric dis- 
tribution around 45°, whereas backward recoiling atoms 
should favor larger angles due to a reduced energy mis- 
match. 

Unfortunately, the analysis of is restricted to short 
times for which the condensate remains practically unde- 
pleted and only first-order side modes are slightly popu- 
lated. As a result there is no evidence that their theory is 
able to explain the characteristic X-shape pattern which 
involves higher-order atomic side modes. On the con- 
trary, there are several different arguments which show 
that the explanation of 0] cannot account for the ob- 
served asymmetry and some of them have already been 
discussed elsewhere Here we would like to give a 

few more arguments which show the failure of the model 
employed in Q . To this end, by analyzing Fig. 1 (A) of 
which depicts a typical experimental outcome, we have 
obtained a rather precise picture of the distribution of 
the atomic side modes which, for the sake of compari- 
son, is presented in Fig. [S] Clearly, the atoms recoil at 
angles much smaller than 45° in the forward direction 
when measured with respect to the BEG center and the 
direction of the applied laser pulse (dashed lines in Fig. 
|3J). The corresponding angle in the backward direction is 
larger than 45°. 

These observations are incompatible with Fig. 3(b) of 
7]. More precisely, this figure shows a very broad an- 
gular distribution for the forward and backward recoil- 



FIG. 5: (Color online) Observed spatial asymmetry in an ac- 
curate schematic representation based on Fig. 1(A) of 0. The 
forward (backward) peaks clearly appear at angles smaller 
(larger) than 45°. The dashed lines indicate the angles at 
45°. 



ing atoms. However, the angular width of the individual 
backwards and forwards travelling atomic wave packets 
observed in the experiments is actually much smaller. 
This means that the curves of Fig. 3(b) should be inter- 
preted as probability distributions for the directions of 
recoiling atoms averaged over many realizations. In this 
case, however. Fig. 3(b) of shows that atoms should 
predominantly recoil at 45° in both forward and back- 
ward directions. Moreover, one should also observe many 
experimental runs where the forward recoiling atoms are 
emitted under a larger angle than the backwards recoil- 
ing ones. None of the experiments, however, reports on 
such types of observations. Furthermore, as depicted in 
our Fig. El experiments show the forward peaks to have 
a much smaller distance from the BEC center than the 
backward peaks. In an explication based on angular dis- 
tributions, however, these distances have to be the same. 

These contradictions between theory and experiment 
clearly indicate that angular effects cannot account for 
the experimentally observed asymmetry and X-shape 
pattern. The main reason for the failure of the model 
employed in 0] is that it does not take spatial effects into 
account which, as we show here, are crucial in producing 
both the asymmetry and the X-shape pattern. Indeed, 
keeping only the two optical endfire modes, one can eas- 
ily verify that in the semiclassical regime the model of 
turns into a spatially-independent mean-field model sim- 
ilar to the ones employed e.g., in and introduced 
at the end of Sec. II. 

As depicted in Fig. Efa), this simpler model cannot 
explain the characteristic X-shape pattern: There is no 
mechanism present to prevent the growth in off-diagonal 
side modes such as (±2,0) and (0, ±2), and they soon 
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FIG. 6: Spatially independent model. Atomic side-mode dis- 
tributions (compare with Fig. (a) Strong-pulse regime: 
tf = 10.4/is and g — 3.0 x 10® s~^; (b) weak-pulse regime: 
tf — 151.5 j-LS and g — 3.5 x 10^ s~^; (c) weak-pulse regime: 
tf = 182 /is and g = 3.5 x 10^ s~\ 

become strongly populated. It is therefore obvious that 
spatial effects play a significant role in the strong-pulse 
regime. 

V. WEAK-PULSE REGIME 

In this section, we will discuss some characteristic dy- 
namic effects arising in the weak-pulse regime of BEC 
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FIG. 7: Spatially dependent model and weak-pulse regime 
(g = 3.0 X 10^ s~^). Spatial distribution of the condensate (full 
curves), side modes (1, 1) (dotted), (1, —1) (dash-dotted), and 
(2,0) (dashed) at times tf = 200 (a); tf = 270 (b); tf = 350 
(c); tf = 420 ^is (d). 



superradiance. Together with Sees. Ill and IV, these 
results illustrate how our model allows to obtain new 
insights into the system behavior. First of all, we de- 
scribe typical stages of the spatially resolved time evolu- 
tion of the matter-wave fields. In particular, we address 
the depletion of the condensate center reported in Ref. 
0|- Secondly, we discuss the time development of the 
optical fields and show that a minimum in the superradi- 
ant light emission does not imply that, at that time, the 
optical fields are also small inside the atomic sample. Fi- 
nally, we show that our model is capable of reproducing 
the characteristic fan patterns for the side-mode distri- 
butions that have been observed experimentally. Addi- 
tionally, we point out similarities and discrepancies in the 
predictions of the spatially dependent and independent 
models. 



A. Time evolution of matter-wave fields 

In this subsection, we examine characteristic features 
of the spatially resolved matter-wave dynamics in the 
weak-pulse regime. We describe the early stages in the 
time evolution of the (1,±1) and (2,0) side modes - 
which are resonantly coupled to the condensate - and 
we discuss the depletion of the condensate center . 

A main difference between the weak- and strong-pulse 
regimes is the fact that in the latter only diagonal side 
modes with \n\ — \m\ are populated significantly. As ex- 
plained in Sec. IV, this is due to the missing spatial over- 
lap between, e.g., the side modes (1, ±1) and the endfire 
mode e± which prohibits the population of the (2, 0) and 
other non-diagonal side modes. Instead, higher-order di- 
agonal side modes are populated, since this is not im- 
peded by spatial effects and the strong fields allow to 
overcome the concomitant detuning barriers. 

A different situation arises in the weak-pulse regime. 
In Figs. we show a typical scenario for the time evo- 
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FIG. 8: Depletion of condensate center in the weak-pulse 
regime {g — 6.0 x 10^ s~^). Spatial distribution of the con- 
densate at times tf = 215 (a) and tf — 350 fis (b). 



lution of the matter- wave fields. Behavior of this kind is 
observed for a large range of values of the external param- 
eters. For very short times [Fig.[7{a)], we have a situation 
similar to the strong-pulse regime, i.e., the first-order side 
modes (1,±1) start to grow at the edges of the conden- 
sate where the optical fields are strongest. The side mode 
(2,0), although resonant with the modes (1,±1), cannot 
yet be populated because the necessary overlap between 
matter-wave and optical fields is missing. However, since 
the applied laser field is weak, higher-order diagonal side 
modes are now populated less efficiently because of the 
detuning barrier. 

Instead, a sort of Rabi oscillation sets in between the 
first-order side modes and the condensate: after the con- 
densate population at some point z has been pumped 
completely to the first-order side mode, it is subsequently 
transferred back. This leads to the appearance of a min- 
imum in the condensate density (marking the point of 
complete transfer) and a concomitant regrowth of the 
population from the edges [Fig.lJfb)]. As to be expected, 
the minima stay near the density maxima of the first- 
order side modes. Once these maxima get close the cen- 
ter of the sample, the overlap between the (1, ±1) side 
modes and the endfire mode e± becomes large, and the 
(2, 0) side mode starts to grow rapidly. Simultaneously, 
the minima in the condensate density merge at the cen- 
ter [Figs.Efc) and [3d)]. This scenario thus shows that 
there is a clear correlation between the onset of popula- 
tion growth in the (2,0) side mode and the depletion of 
the condensate center. We expect that this effect should 
be observable experimentally. One should also note the 
localization of the (2, 0) side mode around the center of 
the system. 

In Ref. 0, depletion of the condensate center was re- 
ported after the system had already undergone a large 
number of superradiant emission cycles and high-order 
side modes were populated. As evidenced in Figs. |S1 our 
simulations are able to reproduce this effect (note that 
g is increased by a factor of 2 compared to Fig. 2} • At 
the times shown in these pictures, side modes of order 
n — 3 and 4 are populated, respectively. We find that 
the depletion of the center is clearly maintained over a 
long period of time. It should be noted, however, that 
the central population increases again briefly after the 
initial depletion described above. The simulations also 
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FIG. 9: (a) Emitted superradiant light intensity as a function 
of time for g — 3.0 x 10^ . Inset: Populations of condensate 
(full curve), side modes (1, ±1) (dashed) and (2,0) (dotted), 
(b) Moduli of total electric field |e+| (full curve) and compo- 
nents le!^'"'! (dashed) and le^'^'l (dotted). 

show that the depletion effect becomes less pronounced 
if the pulse strength is reduced. 

B. Dynamics of optical fields 

Further insight into the characteristics of BEC su- 
perradiance can be obtained from studying the dynam- 
ics of the optical endfire modes. In Fig. Ela) we show 
the emitted light intensity, which is proportional to the 
squared field strengths |e±p at the right and left edges 
of the condensate, respectively. The chosen coupling 
strength g = 3.0 x 10^s~^ is the same as in Fig. [7| 
and the depicted behavior is again representative for the 
weak-pulse regime. For this value of g, only the side 
modes (1,±1) and (2,0) become significantly populated 
in addition to the condensate [see inset in Fig. IHI^a)]. 
To analyze the behavior of the intensity, we use Eq. 
(jl5b|l to decompose, e.g., the endfire-mode field e+ as 
e+(^,^) = E„,™eV"'™)(e,T) with 

e?'"^=-^- f de'e^("-™)^Vn™(?',r)V:+i,™-i(C',r). 

X J — oo 

Each term e^''™'' can be interpreted as the contribution 
to the total field strength e+ arising from transitions be- 
tween the side modes (n, m) and (n -I- 1, to — 1). In Fig. 
Elb), we plot the moduh |e+|, |e^°'"'|, and |e^^'^^|. This 
diagram shows that the first maximum in the radiated in- 
tensity [cf. Fig. Ela)] can be attributed to the transition 
between the condensate and the first-order side modes. 
The second maximum, however, is only partially caused 
by the higher-order transition to the side mode (2, 0), but 
has also a large component due to the continuing Rabi 
oscillations between the BEC and the side modes (1, ±1) 
described in Sec. lV Al 

In Figs. El we show the spatial dependence of the 
optical endfire modes at different times. For some of 




FIG. 10: Spatial distribution of optical endfire modes for g = 
3.0 X 10^ s~^ and times tf = 200 (a); tf = 270 (b); tf = 305 
(c); tf = 420 /is (d). Full curve, |e+j; dotted, |e^'°'|; dashed, 
|e+'"'"''|; dash-dotted, |e_|. 



these times, the corresponding atomic side-mode distri- 
butions are displayed in Figs. |3 Initially, the fields in- 
crease monotonically towards the edges of the condensate 
[Fig. llOf a). compare with Fig. [7fa)]. After some time, 
however, the fields acquire a maximum inside the atomic 
sample [Fig. llOf b). see also Fig.[7fb)]. This is due to the 
fact that near the edges, atoms are transferred back from 
the side modes (1, ±1) into the condensate as mentioned 
in Sec. IV Al This requires the absorption of photons 
from the endfire modes and their subsequent emission 
into the laser field. As a consequence, the endfire-mode 
intensity is reduced. In Fig. [TUT c). the field strength at 
the edges of the condensate, which determines the emit- 
ted superradiant light intensity, has decreased to almost 
zero. However, inside the atomic sample the fields are 
still strong. This means that a small emitted light in- 
tensity does not imply that the fields are also vanishing 
within the sample. Furthermore, from Fig. llOf c) we see 
that the field components e^'^'' and e^'^'' are still strong 
at the condensate edge, i.e., the vanishing of the total 
field strength can be considered as an interference effect. 
For later times, the behavior of the the fields and the 
various components becomes more complex as evidenced 
by Fig. lTUT d) [corresponding to Fig. |3d)]. 



C. Side- mode patterns 

In the experiments on sequential superradiance in the 
weak-pulse regime, the atomic side-mode distributions 
were observed to display a characteristic fan shape [see, 
e.g.. Figs. 1(F), (G) of Ref. ^ and Fig. 1(B) of Ref. '^J. 
As shown by Figs. ^3 and 1121 our theoretical model is 
able to reproduce these patterns. Each of these figures 
displays atomic side mode distributions for a certain cou- 
pling strength and varying pulse durations. Clearly, the 
patterns shown in Figs. |IIi;b)-|nid), IT^b). and [11 d) 
closely resemble the experimentally observed distribu- 
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FIG. 11: Atomic side-mode distributions in the weak-pulse 
regime. The diagrams are for g = 4.5 x 10^ s~^ and pulse 
durations t/ = 140 {a,);tf = 205 (b); t/ = 260 (c); = 335 /is 
(d). Gray level as in Fig. |S1 



(a) 



-10 



12 3 4 

n 




(c) 



3 
2 
1 

^ 
■-1 
-2 
-3 



(d) 



-10 



1 2 3 4 5 

n 



-10 1 2 3 4 5 

n 



arises for short pulse durations before the full fan pat- 
terns develops. In this way, it appears to be intermediate 
between the X and the fan shape. It would be interesting 
to determine whether such patterns can also be observed 
experimentally. So far, however, the side mode distribu- 
tions and their dependence on external parameters such 
as pulse strength and duration have not yet been studied 
systematically in experiments. 

In contrast to the X shape patterns, fan-like side-mode 
patterns can also be obtained with the simpler mean-field 
model using spatial independent amplitudes [see Figs. 
|SJb) and El^c)] . This implies that the role of spatial ef- 
fects for pattern formation is somewhat less significant 
in the weak-pulse regime. Nevertheless, the fan shape is 
not only due to the varying detunings between the dif- 
ferent side modes, but crucially depends on the coupled 
dynamics between the optical and matter-wave fields. 
We conclude this from studies of the mean-field model 
(|16|l with the time-dependent optical endfire modes re- 
placed by constants. In this case, the side-mode dynam- 
ics strongly deviates from the behavior found in the other 
models. 

Although the spatially-independent mean-field model 
can reproduce fan-like patterns, significant differences to 
the spatially dependent model arise when studying very 
weak pulses and long times. These differences are illus- 
trated in Fig. El which shows the time evolution of the 
side mode populations in both models for a small value 
of the coupling constant. In the spatially independent 
model, we see that the population is successively trans- 
ferred in a very regular way from the condensate to the 
side modes (1,±1), (2,0), (3,±1), (4,0), etc. [Figs.UHfa) 
and ll^f b)]. Taking propagation effects into account, 
however, we see that the population remains trapped in 
the side modes (1,±1) and (2,0) for a very long time, 
and the regularity in the time evolution is apparently 
lost completely [see Fig. [CT c)]. 



VI. CONCLUSIONS 



FIG. 12: As Fig. [TTl but for g ^ 7.0 x 10^ s"^ and pulse 
durations tf = 80 (a); tf = 170 (b); tf = 210 (c); tf = 280 /is 
(d). Gray level as in Fig. |S| 

tions mentioned above. Together with Figs. 2(b), (c) of 
Ref. ^3 1 these figures demonstrate that the emergence 
of such patterns is a typical outcome of our simulations 
over a large parameter range in the weak-pulse regime. 

We also point out the similarity of our Fig. I12r c) to 
Fig. 2(c) of Ref. 20]. In particular, both figures show 
population of the side modes(0, ±2) which is not evident 
in Refs. However, we also find a characteristic pat- 

tern which has not been reported experimentally so far 
and in which mainly side modes (n, m) on the "forward 
diagonals", i.e., for |m| = n, n > 0, become populated 
[see Fig. Illf a) and Fig. ^Ja)]. This pattern typically 



We have presented a detailed theoretical analysis of su- 
perradiant Rayleigh scattering from BECs in the frame- 
work of the spatially dependent semiclassical Maxwell- 
Schrodinger equations. Our theoretical approach is 
rather general and equally well applies to many other 
problems involving interaction between light and ultra- 
cold atoms, e.g., matter- wave amplification 0,1121 ■ For 
the problem under consideration, the model allows us 
to reproduce and explain several characteristic features 
observed in the experiments that have not been theoret- 
ically accounted for so far. 

For the strong-pulse regime our studies show that both 
the spatial asymmetry between backwards- and forwards- 
scattered atoms as well as the X-shape side mode pat- 
terns are consequences the spatial inhomogeneity of the 
optical fields. Such effects cannot be explained in the 
framework of spatially independent models. Our results 
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FIG. 13: Weak-pulse regime {g = 2.0 x 10® s"^). Time evolu- 
tion of side-mode populations (a) and distributions (b) for the 
spatially independent model (lltil . (c) Time evolution of side- 
mode populations for the spatially dependent model. Gray 
level as in Fig. |S1 



tween the onset of population growth in the side mode 
(2, 0) and the depletion of the condensate center. In dis- 
cussing the dynamics of the optical fields, we have shown 
that a minimum in the superradiant light emission does 
not necessarily imply that, at the same time, the optical 
fields are also small inside the atomic sample. 

Finally, we have also given an analytical discussion of 
the short-time growth regime of the superradiant process. 
We have shown that the propagation effects turn the ex- 
ponential growth expected from the spatially indepen- 
dent model into a slowed-down subexponential growth. 

Many features and effects that we have discussed 
throughout this work (e.g., the relation between the onset 
of population in the side mode (2, 0) and the depletion of 
the condensate center, or the appearance of intermediate 
side-mode patters) should be experimentally observable. 
We therefore hope that our results will stimulate further 
experimental investigations into the topic. 

There are several directions for further theoretical work 
on the basis of the results presented here. First of all, 
it is of interest to incorporate the effects of the initial 
quantum-mechanical noise, that starts up the process, 
more fully into the theory and to discuss the ensuing 
fluctuations in the macroscopic system dynamics. First 
steps in this direction were reported in [Tl|. Secondly, 
as mentioned above, the present theory can directly be 
applied to study matter-wave amplification. In this way 
it should be possible to obtain an improved control over 
the process. Finally, it would also be important to extend 
the present one-dimensional treatment of the superradi- 
ant dynamics to a full three-dimensional theory in order 
to study transverse propagation effects more closely, such 
as multimode radiation. 



APPENDIX 

We introduce the functions Q(£,,t) = ipi -i{£,,t), 
B(e,T) = ^lia(e,T)e2-, and Eit'r) = 

iLd^'iQi^'j^) + B{C,r)]. Using Voo(e,^) = 
and setting F = F/A, Eqs. (|17|l can be written as 



strongly support an interpretation of the strong-pulse 
regime in terms of atomic scattering from optical fields, 
thus resolving an open controversy as explained in Sec. 
IV. 

In the weak-pulse regime, our model is able to repro- 
duce the experimentally observed fan patterns of the side- 
mode distributions. Moreover, we have also obtained an- 
other kind of pattern which has not been experimentally 
reported so far and which can be regarded as interme- 
diate to the X-shape strong-field pattern. We have pre- 
sented a detailed investigation of the spatially-resolved 
coupled dynamics of matter-wave and optical fields and 
pointed out characteristic effects. In particular, we have 
explained in detail that there is a clear correlation be- 



dQ{^,r) 
dr 

dr 



Q{Lr)+B{tT), (A.la) 
fE{^,r), (A.lb) 
-~TE{^,t) + 2iB{^,t). (A.lc) 



We wish to solve Eqs. (jA.l(l under the boundary condi- 
tions Qi^O) = Vi,-i(?,0), B{^,0) = 0, and E{Q,t) ^ 0. 
Our starting-point is the Laplace transform method ap- 
plied in Ref. |23| to a similar set of equations. We 
thus define the transforms e(s,T) — C{E{^,t)} = 
J^e-^iE{tT)d^, qis,T) = C{QitT)}, and 6(s,t) = 
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C{B{^,t)}. Since E{0,t) = 0, Eqs. l|XT|) transform to 
= q(s,T)+6(s,r), (A.2a) 
= fe(s,r), (A.2b) 

= -fe{s,T) + 2ib{s,T). (A.2c) 



se(s, t) 
dq{s,T) 

dr 
db{s,T) 



dr 



After replacing e(s,T) by (q + b)/s, Eqs. (|A.2bl) and 
(|A.2c|) turn out to have the same formal structure as 
the equations of motion for Ci__i(t) and C-i^i(t) in the 
spatially independent model, and we can immediately 
deduce their solution from Eqs. H18a|l and llSbll by re- 
placing r with 2r/s. This yields 



g(s,r) = 

b{s,T) = 
with Ai(s) = 



2[i-Ai(s)] 

_P^2(s) 



=Ai(s)t 



A2(.s)-- 
s 



Xlis)-- 
s 



q{s,0) 



2[i-Ai(s)] 



i(s) s V J 



(A.3) 
(A.4) 
= i + 



i - ~ 2iT/s and A2(s) 

— 2iT/s. In contrast to the case of Ref. [23l |. 
the inverse transform can no longer be carried out in 
closed form. In view of Eq. (|23|l . we thus specialize 
g(s,0) = l/{s — 2mn/K) which is the Laplace transform 



of the exponential factors. The inversion integral 

-1 nx+ioo 

Q{^,r) = — dse^«g(s,T), (A.5) 

X > 0, is then evaluated by means of the saddle-point 
method (see, e.g., Ref. |2J|). Except for very short 
times, the terms in Eqs. (|A.3|I and (jA.4p involving Ai 
can be neglected, since they only yield an exponentially 
decaying contribution. The relevant saddle point sq for 
the evaluation of Eq. (|A.5|I is then given by the so- 
lution of ^so -|- A2(so)''' = 0. In general, sq has to 
be expressed as the root of a fourth-order polynomial. 
In the weak- and strong-pulse limits of F ^ 1 and 
r ^ 1, however, one can approximate sq ~ vTv and 
So ~ e^*'^/^(f T^/25^)^/^, respectively. Inserting these 
values into the saddle-point formula (and also approxi- 
mating the accompanying slowly varying function appro- 
priately), we obtain relations (fH^ . (|?T|l . and f^ . We 
have verified that the saddle-point method without these 
further approximations as well as - within their respec- 
tive ranges of validity - the simplified formulas H19|) , H21|l , 
and H24|) give very good approximations to the numerical 
solution of Eqs. (|17|l . The result H19|) is consistent with 
the fact the solution of Eqs. H17|l without the backward 
mode 'ip-i^i{S,,T) is given by ?Ai,_i($,r) = V'o/o(2V Tv|) 
with the modified Bessel function Io{x) ~ exp(a;)/'\/27ra, 
a; S> 1 (see also 0)- 
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